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Following the earlier approach [J| the explicit expression for the two time free energy distribution 
function in one-dimensional random directed polymers is derived in terms of the Bethe ansatz replica 
technique. It is show that such type of the distribution function can be represented in the form of 
a generalized ’’two-dimensional” Fredholm determinant. 
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I. INTRODUCTION 

The problem of directed polymers in a quenched random potential or equivalent problem of the KPZ-equation Q 
describing the growth in time of an interface in the presence of noise have been the subject of intense investigations 
during past three decades Id hla. A few years ago the exact solution for the free energy probability distribution 
function (PDF) has been found [Ibl-I^. It was show that this PDF is given by the Tracy-Widom (TW) distribution 
[^ . The two-point free energy distribution function which describes joint statistics of the free energies of the directed 
polymers coming to two different endpoints has been derived in [30l432l| . 

In all these studies, however, the problems were considered in the so called ’’one-time” situation. For the first 
time the joint probability distribution function for the free energies of two directed polymers with fixed boundary 
conditions at two different times has been studied in the paper [l| (in terms of non-rigorous replica Bethe ansatz 
approach) and in the paper Q (mathematically rigorous derivation). Unfortunately, the results obtained in these 
papers are somewhat inconclusive: the final formulas are expressed in terms rather complicated mathematical object 
whose analytic properties are not known. Moreover, for the moment it is even not clear whether these two results 
obtained it terms of two different approaches do coincide. 

In this paper I am going to derive the joint probability distribution function for somewhat different two time free 
energy object. It turns out that compared to the previous calculations ’’more symmetric” structure of this 

object (see below) makes the corresponding calculations much more simple and the final result can be expressed in 
the form of a compact formula which could be conditionally called a ’’two-dimensional” Fredholm determinant with 
explicit expression for its kernel, eqs. (EHi-dMl). 

The model under consideration is defined in terms of an elastic string (/)(t) directed along the r-axes within an 
interval [0, t] which passes through a random medium described by a random potential V ((/>, r). The energy of a given 
polymer’s trajectory (fir) is 


where the disorder potential r] is Gaussian distributed with a zero mean = 0 and the ^-correlations 

V{(j),T)V{(j)' ,t') = u6{t — T')S{(j) — 4>')- The parameter u describes the strength of the disorder. For the fixed 
boundary conditions, (fiO) = 0, (j)(t) = x, the partition function of this model is 

Zt{x) = / T>(/)(t) = exp{-l3Ft{x)) (2) 

J</i(0)=0 

where /? is the inverse temperature and Ft{x) is the free energy. In the limit t ^ oo the free energy scales as 
PFt{x) = fffot + fix'^/2t + Xtf{x), where /o is the selfaveraging free energy density, At = oc and /(x) 

is a random quantity described by the Tracy-Widom distribution [2ll - l24| . As the first two trivial terms of this free 
energy can be easily eliminated by simple redefinition of the partition function, they will be omitted in the further 
calculations, in other words, Zt{x) = exp(—At/(x)). 

Let us consider two partition functions: the first one, Zt{0), is defined in eq.(I2]) with the zero boundary conditions, 
(^(0) = = 0, while the second one, ^t.At(O), in addition to the above zero boundary conditions contains an 
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additional constrain (j){t — At) = 0. In other words, in the second case the polymer trajectory before coming to zero 
at time t is forced cross the zero at some intermediate time (< — At). 



FIG. 1: Schematic representation of two directed polymer paths: (1) with <^(0) = (j>{t) = 0, and (2) with (^(0) = (j>{t — At) = 
4>{t) = 0 


By definition, 


/ + 00 

dx Zt-At{x) Z^tix) = exp(-At/) 

-OO 


Zt(0) = 


and 


ZtMO) = Zt-At{0) Z*^M = exp(-At/) 


( 3 ) 


( 4 ) 


where Z^^(x) is the partition function of the directed polymer system in which time goes backwards, from t to (t — At). 
For technical reasons (for proper regularization of the integration over cc at ± infinities) it is convenient to split the 
partition function Zt{0) into two parts, the ’’left” and the ’’right” ones: 

/ O ^+oo 

dx Zt-At{x) Z’^^ix) + / dy Zt-At{y)Z*^^{y) (5) 

-OO J 0 

In terms of the free energy definitions ([3])-(j4]) one would like to compute the joint probability distribution function 

A) = lim Prob[/ > /i; / > / 2 ] (6) 

where it is assumed that in the limit t ^ oo the parameter A = At/t remains finite. In terms of the above partition 
functions ©-O this quantity can be defined as follows: 


kF(/i,/ 2 ,A)= hm ^ ^ 

f .—^on f ^ ^ 


(- 1 ) 


N+K 


t—¥OC 

N^OK^O 


■ exp(AtiV/i + \tKf 2 ) Zf (0) [Zt_At{0) Z- 


1 K 


At'. 


( 7 ) 


where (...) denotes the averaging over the random potential. Taking into account eq.(IS]) one gets 


IT(/i,/ 2 ,A) = lim ^ 


(- 1 ) 


L+R+K 


t^oo ^ L\R\K\ 

L.R.K=0 


■ exp 


Xt{L + R)h + \tKf2 


(* + 00 


dxi...dxL / dyR...dyi x 


' L R \ / L R 

[n^t-At(a:a)]z,^^,(0)[n^t-At(2/6)] [n^At(^a)]z^,^(0)[n^lt(y6) 


6=1 


6=1 


( 8 ) 


Introducing: 


'S{xi, ...,XN;t) = Zt{xi) Zt{x2) ■■■ Zt{XN) 


( 9 ) 
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one can easily show that 4'(x;t) is the wave function of iV-particle boson system with attractive ^-interaction: 


N N 

I3dt'b{x;t) = - + - n'^Sixa - Xt) «'(x;t) 


( 10 ) 


a^b 


(where k = P^u) with the initial condition 'I'(x; 0) = n^]^(5(a;a)- The wave function ^'(x; t) of this quantum problem 
can be represented in terms of the linear combination of the corresponding eigenfunctions of eq. OT. A generic 
eigenstate of such system is characterized by N momenta {Qa} (a = which split into M (1 < M < N) 

’’clusters” each described by continuous real momenta qa {a = and characterized by Ua discrete imaginary 

’’components” (for details see 


IK , 


Qa ^ Qr = 9a - y(?T-a + 1 - 2r) ; (r = l,...,n„; a = l,...,M) 


with the global constraint Explicitly, 


N 


^qw=^ n 


V l<a<b 


. sgn(a;a - xQ 

^"Qv.-Qn 


exp 


N 

^iy^^Qv^Xg 


( 11 ) 


( 12 ) 


where the vector Q denotes the set of all N momenta eq. ED and the summation goes over A^! permutations 7^ of 
momenta Qa, over N particles Xa- In terms of the above eigenfunctions the solution of eq. (nnD can be expressed as 
follows: 


'I^(x;t) = ^ l^^(Q)l' 'I'q(x)^q( 0) exp(-tA(Q)) 


(13) 


where the symbol f Vq denotes the integration over M continuous parameters {gi, the summations over M 

integer parameters {ni,..., um} as well as summation over M = \,N. |C'(Q)p is the normalization factor, 


|C(Q)|2 = 




M I ^ I ^ 

Qoc- q&- - n0)\ 


n 


na=l('«^a) a^P 19a “ 9/3 - "-jijla + np)\ 

1 


det 


^KUa — iqa + ^Knp + iqp 


(14) 


and A(Q) is the eigenvalue (energy) of the eigenstate 'I'q(x) 


N 


M 




M 




(15) 


06—1 

In terms of the wave functions ©-(USD the probability distribution function (0 can be expressed as follows: 


IT(/i,/ 2 ,A) = hm ^ 

T .—Vrvi « * 


t^oc ^ L\K\R\ 

L,K,R^0 

r-\-oo 


exp 


Xt{L + R)fi+XtKf2 


(16) 


/ U r+oo 

dxi...dxL J {t-At)) ---,2/1; At) 


K 


K 


where the second (conjugate) wave function represent the ’’backward” propagation from the time moment t to the 
previous time moment t — At. Schematically the above expression is represented in Figure I. The above expression 
is quite similar to eq.(23) in [l| although here its structure is ’’more symmetric” which essentially simplify further 
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FIG. 2: Schematic representation of the directed polymer paths corresponding to eq. (Hi 


calculations. Repeating the same steps as in the paper [I|, in the limit t —>■ oo one eventually gets (sf eq.(65) in Q): 


Wih,f2,A) = 1 +J 2 


M-. 


(-1)^ 

- (M!)2 


M 

n 


+ 00 


dqadpo 

(27r)2 


/ +00 

dy Ai{y + (1 - A)ql + Apl - /i) x 

-oo 


X [ dziadz^adz^J --'5(^1 a)^(2:2a )^(^3a “ Pa , , ^2a , X 

X exp{(zi„+22a + ^3a)2/ + 23a(/2-/l)} 


X det 


1 


•21a + •22a + •23a “ *9a + Zla' + Z2a' + 23a' + *9a' 

1 


X det 

2la + 22a + 23a + Wa + 2i„/ + Z2a' + 23a' “ Wo 
where the integration contour C is shown in figure 2, and 

2l 




;' = 1.M 


j(q - P,Zl,Z2,Z3^ ^ + 


1 + 


22 


Z2 + Z3 + ^i(q - p)(“) / I " ' Z1 + Z3 - ^i(q - 


(17) 


(18) 


Here (q — p)^^^ = (p — q) A ie where the parameter e has to be set to zero at the end. Determinants in eq. (HB can be 
represented in the following integral form: 


det 


1 


2a - iqa + zp + iqp 


M 


= X! fl “ *'2'“ + + '^dVc) 

a,/3— 'P^Sm 


a=l 


M 


n 

a—l 

M 

n 

a=l 


dUn 


dUa 


(-l)M exp|-(zQ, + Zp^)ua + ^(ga - g'P„)Ma| 

’ESm a=l 

M 

(-l)^^-l exp|-0„(Ma + + *ga(Ma - 


J -PGSm 


a=l 


(19) 
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FIG. 3: Contour of integration C in eq. (HU. 


Substituting this representation into eq. (HU and performing integrations over zi , Z 2 and Z 3 one obtains the following 
result 


B3A./.A) = 1 : (-i)W+Mn 


M 


- (M!) 


V,V'&Sm 




dUadVn 


M 

n 


+00 


dqadpa 

(27r)2 


f dy AiU + (1 - A)ql + Ap^ - fi + Ua + Va + u-p^ + 

J —00 


X exp|+jga(ua - Up^) - ipa{va “ S[{qa - Pa)', V, (/2 “ /l)] 


( 20 ) 


where 


S[{q-p)', y, (/ 2 -/ 1 )] = 47r(5(t/)(5(q-p) 6»(/i -/a) + 


+ 


0 if 2 -fi + y) 0 if 2 -fi-y) - 4%) 


sin[4(g -p )(/2 - /i)] 

(q-p) 


9{f2 - fi) (21) 


According to the above equation in the cases: (a) /2 < /i and (b) /2 > /i, one finds two essentially different results 
for the probability distribution function W{fi,f 2 ,A). 

(a) /2 < /i. In this case the distribution function W{fi,f 2 , A) turns out to be independent of /2 and A: 




M 


- (M !)2 


v,v'eSM 


dUadVn 


M 

n 


/ — Ai^g^ - /l + {Ua +VV'J + (Va + W-p^)) expUqo,[{Ua + Vp>J - {Va +M'P„)]} 

y —00 TT \ / r j 


= 1 + -—— 

M=1 ^ ' V,V'&Sm 


Y, (-i)Md^ln 


OL — 1 


dUadVa 


M 

n 

a—l 


22/3 Ai[2i/3(y^ + _ fj2))] Ai[2^/^{up^ +va- fi/2))] 


( 22 ) 
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Redefining ^ 2 Va ^ 2 and taking into account that 


M 


M 


JJ Ai uv^ -\-Va- = JJ Ai Ua + v^-1 - /l/2^/^ 


a=l 


Ct — 1 


we obtain 


«'(/■) = !+ElMTirll 

= 1 ^ Q!=l 


M 


dva 


(-ijH+Mfj 


V,V'^Sm 


OL — 1 


K{v-p-i-hf2^/^-, vv^-h/2^1^) 


where 


pOO 

K[v\ v') = du Ai(u + ?;) Ai(u + ?;') 

Jo 

is the Airy kernel. Redefining the permutations, V + V' —>• V, we eventually get 


(23) 


(24) 


(25) 


»'(/■) "I+eSswII 

M=1 a=l 


dVn 


M 

E (-i)W n 


VgSm 


np„-/i/22/3) 


= F2{-fil2^/^) (26) 


Thus we have got the usual Tracy-Widom distribution for the free energy /i of the directed polymer with the zero 
boundary conditions, as it should be (it is evident that in the limit t —^ oo with the probability one the free energy 
(/ 2 ) of the polymer which is forced to pass par the zero at some intermediate time (t — At) is larger than the free 
energy (/i) of the polymer which is free of this condition). 

(b) /2 > /i- In this case the situation becomes more tricky. According to eas. d^ and (1^ we get: 


W{fuf2,A) 


/2>/l 


/ I'.M M 

= 1 + -—— TT 

^ (M !)2 li 

M=1 ' ' a=l 


dUndVn 


Y (-i)H+Mn 


v,v'eSM 


where 


d^(Ua (^-Pa a ^ 

(27) 


d^(u,v)‘, {u\v')(^flT f2: A) J' 


(/2—/l) 


dy 


-(/2-/l) 


+ 00 


dqdp 


Ai [y - /i + (1 - A)g^ + +U + V + u' + v'] 


X exp{ig(u — m') — *p(z; —'c')} — 

dqdp sin(i(g-p)(/2-/i)) 


-4 


+ 00 


( 2^)2 


q-p 


Ai [—/i + (1 — A)g^ + Ap^ + uV + u' + v'^ x 

(28) 


X eyip{^iq{u — u') — ip{v — v')^ 

The structure in eg. (1271) could be conditionally represented as a ’’two-dimensional” Fredholm determinant: 


W{h,h,A) 


/2>/l 


= det 1 - K\ 


(29) 


where K = (/i, / 2 , A) is the integral operator function of three parameters /i, /2 and A on the two-dimensional 

space ^ = (u, v) with (rt, v) € [0, oo) and with the kernel given in ea. (l28t . 

Eas. (l77l) - (l^^ constitute the central result of this work for the two-time free energy probability distribution function 
of two directed polymers with the two types of the zero boundary conditions: ^(0) = = 0 and ^(0) = (j){t — At) — 

(j)(t) = 0 in the thermodynamic limit t ^ oo such the parameter A = At/t remains finite. Analytic properties of the 
mathematical object defined in Eas. (l27l) - (l28l) remains to be investigated. 
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